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Abstract 

We present measurements of the flavor-changing neutral current decays B — > Kl~^£~ and B 
K*i~^i~, where is either an e~^e~ or fj,~^i-i~ pair. The data sample comprises 229 x 10^ 

T(45) BB decays collected with the BABAR detector at the PEP-II e^e~ storage ring. We 
measure the branching fractions 

13{B Kl+r) = (0.34 ± 0.07 ± 0.03) x lO"*^ 

BiB K*£+i-) = (0.78l°:{? ± 0.12) x 10~^ 

the direct CP asymmetries of these decays, and the relative abundances of decays to electrons and 
muons. 
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1 Introduction 



The decays B K£^i~ and B K*£~^£~ , where i^£~ are the charged lepton pairs e~^e~ or n~^fJ-~ 
and K* is the K*{892) meson, result from b ^ s flavor-changing neutral currents (FCNC). In 
the Standard Model (SM) of electroweak interactions, such 6 — > s processes are forbidden in tree- 
level Feynman diagrams; they are allowed at lowest order through one-loop diagrams involving the 
emission and re-absorption of W bosons. Because the lowest order SM diagrams are loops of weakly 
interacting particles with virtual energies comparable to the electroweak scale, new flavor-changing 
interactions at the electroweak scale can introduce loop diagrams with comparable amplitudes. 
The SM predictions of the rates and kinematic distributions of FCNC decays can be significantly 
modified by a broad class of new physics models, such as a charged Higgs boson [1], topcolor [1], 
weak-scale supersymmetry [2,3], fourth-generation fermions [4], or leptoquarks [5]. 

In the SM, three amplitudes contribute at lowest order to the b si~^i~ process: a photon 
penguin, a Z penguin, and a W~^W~ box diagram (Figure P). The magnitude of the photon 
penguin amplitude is well known experimentally from measurements of the rate of the FCNC 
decay 6 — > 57 [6] and agrees well with SM predictions [7]. The latter two amplitudes are not well 
known and thus studies of 6 — > s£~^£~ provide new information on FCNC processes. The SM decay 
rate of 6 — > s£~^£~ is suppressed relative to other b decays, resulting in a predicted total branching 
fraction of (4.2 it 0.7) x 10^® [3], in agreement with experiment [8]. 

The most abundant exclusive decays associated with the b — > si'^i^ transition, B K£^i~' and 
B — > K*£~^£~ , are predicted to have branching fractions of 0.4 x 10~^ for B — > K£~^£~ and about 
three times that for B K*i~^£~ [3,9-13], with a theoretical uncertainty of 30%. The theoretical 
uncertainty is predominantly due to the uncertainty in the prediction of semileptonic form factors, 
which model the rate that a 6 ^ s FCNC in a -B decay results in a sing le K^*^ meson. The 
partial widths of S — > Ke^e~ and B — > Kfi^^~ are expected to be identical, because of identical 
electroweak couplings of electrons and muons. The branching fractions of both B K*e~^e~ 
and B K*fj,^n~, however, receive a contribution from a pole in the photon penguin amplitude 
at = m'j+^- ~ 0, and the enhancement in the electron mode is significantly larger due to its 
lower threshold. This phase space difference in the pole contribution is expected to reduce 
the ratio T(B — > K* 11'^ jx') /T{B — > K*e^e^) from unity to 0.752 [3]. Previous measurements 




q q q q 

Figure 1: Examples of Standard Model diagrams for the decays B K^*^£~^i~ . For the photon or 
Z penguin diagrams on the left, boson emission can occur on any of the b, t, s, or W lines. 
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of the exclusive decays are consistent with predictions [14-16]. In the absence of new physics 
contributions, improved precision in the exclusive branching fractions will improve experimental 
constraints B ^ K^*^ form factors. 

More precise SM tests can be obtained from rate asymmetries and kinematic distributions of 
the exclusive decay products. The direct CP asymmetries 

_r{B ^ K^*k+i-) -r{B ^ K(*k+i-) 

~ T{B K(*)£+i-) + r{B Ki*)£+i~) 

for these decays are expected to be very small in the SM, much less than 1% [17], whereas new 
physics at the electroweak scale could enhance Acp to values of order one [18]. If one neglects the 
pole region (g^ < 0.1 GeV^/c^) of B ^ K*e~^e~, in the SM the ratios 

Rk = B{B Kfi+^x-)/B{B Ke+e-), 

Rk* = B{B K*fi^fi-)/B{B K*e+e-) 

are expected to be unity with high precision. However, this ratio could be enhanced by corrections 
of order 10% due to the presence of a supersymmetric neutral Higgs boson with large tan j3 (ratio of 
vacuum expectation values of the two Higgs doublets) [19]. The Feynman diagram of this process is 
shown in Figure EJ A large tan (5 would enhance the Higgs / squark coupling, and the Higgs decays 
to muons will be enhanced relative to decays to electrons because of the large ratio of Yukawa 
couplings m^/rrig. A measurement of the relative abundance of electrons and muons in exclusive 
decays is therefore a probe of scalar penguin processes and complements the limits obtained from 
searches for the rare decay Bg [20]. With sufficiently large samples of i3 ^ K*i~^£~ 

events, angular asymmetries in the four-particle final state can also accurately gauge the relative 
phase and magnitude of the three contributing FCNC amplitudes [2,3, 17]. 

2 Detector and Datasets 

We analyze data collected with the BABAR detector at the PEP-II storage ring at the Stanford 
Linear Accelerator Center. The data sample comprises 208.0 fb"^ recorded on the T{4S) resonance, 




X 



Figure 2: Feynman diagram of a Higgs penguin process which would enhance b sfi^fi relative 
to 6 — > se~^e~ . 



9 



yielding (229.0 it 2.5) x 10^ BB decays, and an off-resonance sample of 22.1 fb~^ used to study 
continuum background. 

The BABAR detector is described in detail elsewhere [21]. The most important capabilities of the 
detector for this study are charged-particlc tracking and momentum measurement, charged n/K 
separation, and Icpton identification. Charged particle tracking is provided by a five-layer silicon 
vertex tracker (SVT) and a 40-layer drift chamber (DCH). The DIRC, a Cherenkov ring-imaging 
particle- identification system, is used (along with dE/dx measured in the trackers) to separate 
charged kaons and pions. Electrons arc identified using an electromagnetic calorimeter (EMC), 
which comprises 6580 thallium-doped Csl crystals. These systems are mounted inside a 1.5 T 
solenoidal superconducting magnet. Muons are identified in an instrumented flux return (IFR), in 
which resistive plate chambers are interleaved with the iron plates of the magnet flux return. 

Simulated samples of signal B decays, charmonium B decays, generic BB decays, and continuum 
e"'"e~ — > qq (for q = u, d, s, or c) events are used to compute selection efficiencies, optimize 
event selection, and estimate certain backgrounds, as described below. The simulation is based on 
GEANT4 [22] detector emulation software. The model for simulating signal B decays is a 6 — s£'^i~ 
matrix element calculation, which includes 0{as) and 0{AQQY^/mb) corrections [3], convolved 

with B —>■ K^*^ form factors predicted by light-cone QCD sum rules [10]. 

3 Event Selection 

We select events that include two oppositely charged lepton candidates {e'^e~, ji^ji'), a kaon 
candidate (either or K^), and, for the B K*i~^£~ modes, a vr^ candidate that, when 
combined with a kaon candidate, forms a K* candidate. Electron (muon) candidates are identified 
by a likelihood (neural-net) based algorithm, and are required to have a minimum momentum 
p > 0.3GeV/c {p > 0.7GeV/c) in the laboratory frame. 

Bremsstrahlung photons from electrons are recovered by combining an electron candidate with 
up to one photon with > 30 MeV. Recovered photons are restricted to an angular region in 
the laboratory frame of {9^, 4>^) = {6e ± 35 mrad, (pe ± 50 mrad) around the initial electron direction 
{9e,<pe)- Photon conversions and 7r° Dalitz decays are removed by vetoing all e+e~ pairs with 
invariant mass less than 0.03GeV/c^, except in B ^ K*e^e~ modes, where we preserve acceptance 
at low invariant masses by retaining pairs that intersect inside the beam pipe. 

Charged kaon candidates are tracks with dE/dx and DIRC Cherenkov angle consistent with 
the angle expected for a kaon. tt^ candidates are tracks that do not satisfy the selection. K'g 
candidates are reconstructed from two oppositely charged tracks with an invariant mass (computed 
assuming they are vr+vr") consistent with the mass and a common vertex displaced from the 
average interaction point by at least 1 mm. 

True B signal decays produce narrow peaks in the distributions of two kinematic variables, which 
can be fitted to extract the signal and background yields. For a candidate system of B daughter 
particles with total momentum pB in the laboratory frame and energy E"^ in the ^(45*) center- 

of-mass (CM) frame, wc define mss = \l (•''/2 + c^po ■ Pb)^/-E'o ~ ^^Ps <^iid /S.E = E*^ — \/s/2, 
where Eq and po are the energy and momentum of the T(4S') in the laboratory frame, and y/s 
is the total CM energy of the e+e~ beams. For signal events, the mss distribution peaks at the 
B meson mass with resolution a ~ 2.5 MeV/c^, and the distribution peaks near zero, with 
a typical width o" ~ 20 MeV. In i? ^ Ki'^i~ channels, we perform a two-dimensional unbinned 
maximum- likelihood fit to the distribution of ities and AE in the region mES > 5.2 GeV/c^ and 
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|A£^| < 0.25 GeV. In B ^ K*i~^£~ decays, we perform a three-dimensional fit to rriES, in the 
same regions as for B — >^ Kl'^l~ , and in addition we include in the fit the kaon-pion invariant mass 
for the region 0.7 < rriKir < 1-1 GeV/c^. 

Backgrounds arise from four main sources: (1) random combinations of particles from qq events 
produced in the continuum, (2) random combinations of particles from T(4S') BB decays, (3) B 
decays to s£'^£~ final states other than the signal mode ( "crossfeed" ) and (4) B decays to topologies 
similar to the signal modes. The first two ("combinatorial") backgrounds typically arise from pairs 
of scmilcptonic decays of D or B mesons and produce distributions in mss and AE which are 
broadly distributed compared to the signal. The third source has niES similar to signal, but the 
peak of the AE distribution is significantly offset from the signal due to the addition of a random 
particle ("feed-up") or omission of one of the B daughters ("feed-down"). The last source arises 
from modes such as B ^ J/ijjK^*^ (with J /tj) t^l~) or B ^ K^*^TnT (with pions misidcntificd as 
muons), which have shapes similar to the signal. All selection criteria are optimized with simulated 
data or with data samples outside the region of the maximum-likelihood fit. 



3.1 Combinatorial backgrounds 

We suppress combinatorial background from continuum processes using a Fisher discriminant [23] , 
which is a linear combination of variables with coefficients optimized to distinguish between signal 
and background. The variables used in the Fisher discriminant are the following kinematic quanti- 
ties computed in the CM frame: (1) the ratio of second- to zeroth-order Fox- Wolfram moments [24] 
for the event, computed using all charged tracks and neutral energy clusters; (2) the angle between 
the thrust axis of the B candidate and that of the remaining particles in the event; (3) the pro- 
duction angle 6b oi the B candidate with respect to the beam axis; and (4) the masses of Ki pairs 
with the same charge correlation as a semileptonic D decay. The first three variables exploit the 
differences in event shapes between the jet-like topology of light quark pair production and the 
spherical shape of T{4S) — > BB production. The fourth variable discriminates between D meson 
decays, which have Ki mass distributed below the D mass, and signal decays, which have a broad 
distribution in Ki mass. 

We suppress combinatorial backgrounds from BB events using a likelihood function constructed 
from (1) the missing energy of the event, computed from all charged tracks and neutral energy 
clusters; (2) the vertex fit probability of all tracks from the B candidate; (3) the vertex fit probability 
of the two leptons; and (4) the angle 9b- Missing energy provides the strongest suppression of 
combinatorial BB background events, which typically contain neutrinos from two semileptonic B 
decays. 

The parameters of the Fisher discriminant and the likelihood function are determined separately 
for each of the eight signal decay modes. The selection criteria for the background suppression vari- 
ables are optimized simultaneously, and are chosen to minimize signal yield statistical uncertainties 
in each mode. The efficiencies of the Fisher and likelihood requirements are validated by comparing 
the efficiencies in data and in simulation using the B — > J/tpK^*^ control sample. 



3.2 Peaking backgrounds 

The largest backgrounds that peak in uies and AE are B decays to charmonium: B J/ipK^*') 
(with J/V' ^ i+i-) and B ^{2S)K^*^ (with ^(25) ^ i+i'). We exclude dilepton pairs 
consistent with the J/tp mass (2.90 < rrig+e- < 3.20 GeV/c^ and 3.00 < m^+^- < 3.20 GeV/c^) or 
with the ip{2S) mass (3.60 < m^+£- < 3.75 GeV/c^). This veto is applied to m^+^- both with and 
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without the inclusion of bremsstrahlung photon recovery. When a lepton radiates or is mismeasured, 
m£+£- can shift away from the charmonium mass, while AE shifts in a correlated manner. The veto 
region is extended in the {m£+£- , AE) plane to account for this correlation (Figure EJ, removing 
nearly all charmonium events and simplifying the description of the background in the fit. Because 
the charmonium events removed by these vetoes are so similar to signal events, these modes provide 
large control samples (about 13700 events of S — > J/tpK^*^ and 1000 events of -B ^ ■0(25')i^^*^ 
in all) for studying signal shapes, selection efficiencies, and systematic uncertainties. After the 
vetoes on B J/tpK^*^ and B '4){2S)K'^*^ decays, the remaining peaking background from 
these processes is estimated from simulation to be 0.0-1.6 events, depending on the decay mode. 

In the muon modes, the pion misidentification rate is significant (~ 2%), leading to additional 
peaking backgrounds from the decay B~ D^tt^ with K^tt^ or K*^tt^, or from 

B^ D+TT- with D+ K*^TT+ (and their charge conjugates [25]). These events are suppressed 
by vetoing events where the X*-*-'// mass is consistent with a hadronic D decay. The remaining 
background from the charmless decays B K^*>TnT, B K^*^K7r, and B K^*^KK is estimated 
from data. We select control samples of i? ^ K^*^hfi events with the same requirements as signal 
events, except that muon particle identification is no longer required for the hadron candidate h 
and hadron identification requirements for pions and kaons are used instead. This results in a 
sample of predominantly hadronic B decays. Each event is given a weight corresponding to the 
muon misidentification rate for the hadron divided by its hadron identification efficiency, and the 
number of peaking background events from hadronic B decays is extracted from the weighted mES 
distribution through a maximum-likelihood fit similar to that used to extract a signal. The 
control sample calculations result in an estimate of 0.4-2.3 background events per decay channel 
for B — > K^*^ modes. 

Finally, there is a peaking contribution to the electron modes from the rare decays B K*^ 
(with photon conversion in the detector), B K^*^-k^, and B K^*^r] (with a vr*^ or rj Dalitz 
decay to e~^e~'y). The sum of these backgrounds is estimated from simulation to be 0.0-1.4 per 
decay channel for the B — > K^*^e~^e~ modes. 

The number of peaking background events from all sources is shown in Table ^ for the in- 
dividual decay modes. The peaking backgrounds in the modes with electrons are dominated by 
processes with real electrons, and the uncertainties are dominated by simulation statistics. The 
peaking backgrounds for modes with muons are dominated by hadrons misidentified as muons; the 
dominant uncertainty here is systematic and originates from the unknown K/tt composition of the 
contributing hadrons. 
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Figure 3: Charmonium veto regions in the — K'^e'^e~ (above) and S+ — ji^ ji^ (below) 
channels. The points are simulated J/V' and i^i^S) events, with abundance equal to the mean 
number expected in 208 fb~^. The projections onto m£+£- and A£' are shown above and at 
left, which indicate the high density of points at {mi+i- , A.E) = (m^,0.0). The vertical band 
corresponds to events where the J/ip {ip{2S)) and K'^ come from different B decays; the diagonal 
band corresponds to events with mismeasured leptons. 
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Table 1: Mean expected peaking backgrounds in 208 fb ^, for the individual K^*^^'^^ decay modes 
after applying all selection requirements. 

Decay mode Events 

B+ K+e+e- 0.7 ±0.2 

B+^K+n+fi- 2.3 ±0.5 

K^e+e- 0.01 ± 0.01 

B^ K^^+^j,- 0.4 ±0.1 

^0 _^ K*^e+e- 3.0 ±0.6 

^ K*^^+n- 1.4 ±0.8 

B+-^K*+e+e- 0.9 ± 0.2 

B+^K*+n+n- 0.6 ±0.3 
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4 Fits 



For B Ki^£~ , a two-dimensional fit to tit-es and AE is performed. For B K*l~^£~, the mass 
of the K* is added as a third fit variable. The signal shapes are parameterized with separate Crystal 
Ball functions [26] for rriES and AE. Both the uies and AE shape include a radiative tail, which 
accounts for the effects of bremsstrahlung of the electrons in the BABAR detector. The mES shape 
parameters are additionally assumed to have AE dependence cq + C2{AE)'^; the variation of the 
rriES width due to the quadratic term is typically a few percent of cq. All signal shape parameters 
are fixed from the signal simulation, except for the mean and width parameters in m-Es and AE, 
which are fixed to values from charmonium data control samples (for the rriES width, cq is fixed 
from charmonium data and C2 is fixed from signal simulation). In the B K*i~^i~ channels, the 
K* is fitted with a relativistic Breit-Wigner line shape. Adding the mass of the K* to the likelihood 
fit increases the precision of the B — > K*i~^£~ branching fraction measurement by approximately 
10%. 

The background is modeled as the sum of three or four terms: (1) a combinatorial background 
shape with floating normalization, written as the product of an ARGUS function [27] in rwES) a 
linear term in AE, and the product of ^Jrax-K — 'mx — rrin and a quadratic function of mKn for 
the K* modes; (2) a peaking background contribution, with the same shape as the signal, but 
with normalization fixed to estimates of the mean peaking backgrounds (see Table and (3) 
terms with floating normalization to describe (a) background in B ^ Kl^l~ (B K*i~^i~) from 
B K*e+£- {B K*7r£+e-) events with a lost pion, and (b) background in B ^ K*£+£- 
from B K£^i~ events with a randomly added pion. In the K* modes, we allow an additional 
background (4) that uses our combinatorial shape in rriEs and AE, but peaks in rriKn at the 
K* mass. The yield of this term is fixed to (5 it 5)% of the total combinatorial background, as 
determined from simulation. Because the normalizations for terms (1) and (3) are floating, as are 
the combinatorial background shape parameters, much of the uncertainty in the background is 
propagated into the statistical uncertainty on the signal yield obtained from the fit. 

The direct CP asymmetry Acp is also extracted from the fit to the modes B~^ — > K~^i'^£~ and 
B K*£^£~ , where the h flavor of signal candidates can be inferred directly from the charges of the 
final state K^*'> hadrons. It is not possible at this time to measure Acp in the mode B^ — > K^£^£~ , 
as the signal statistics are small and the h flavor can only be inferred indirectly from properties of 
the other B meson. The CP asymmetry of the combinatorial background is allowed to float in the 
fit, while the asymmetries of the peaking background and crossfeed background are fixed to and 
varied from -1 to 1 to evaluate the systematic uncertainty associated with these components. 

5 Systematic Uncertainties 

Table 121 lists the relative systematic uncertainties on the efficiency for each mode. The sources of 
uncertainty considered are: charged-particle tracking (0.8% per lepton, 1.4% per charged hadron), 
charged-particle identification (0.5% per electron pair, 1.3% per muon pair, 0.2% per pion, 0.6% per 
kaon), the continuum suppression cut (0.3%-2.2%, depending on the mode), the BB suppression 
cut (0.6%-2.1%), selection (0.9%), signal simulation statistics (0.4%-0.7%), and the number of 
BB events (1.1%). The uncertainty in the signal efficiency due to model dependence of form factors 
is evaluated for each mode to be the full range of variation from a set of models. The models con- 
sidered are based on QCD sum rules [9], light-cone QCD sum rules [10], and lattice QCD [11]. The 
model dependence enters through the variation in distributions; since the selection efficiency is 
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not highly sensitive to this distribution the efficiency varies by only 4%-7%. For branching fraction 
measurements which combine modes, the systematic uncertainty is an appropriately weighted sum 
of correlated and uncorrelated sources from the contributing modes. The total systematic uncer- 
tainty in the signal efficiency introduces a systematic uncertainty A;Bg|j in the measured branching 
fraction. 



Table 2: Systematic uncertainties of selection efficiency (in %) considered for K^*^i'^£ decays. 



Source 


















TrkefT. {e, fi) 


±1.6 


±1.6 


±1.6 


±1.6 


±1.6 


±1.6 


±1.6 


±1.6 


Electron ID 


±0.5 




±0.5 




±0.5 




±0.5 




Muon ID 




±1.3 




±1.3 




±1.3 




±1.3 


Kaon ID 


±0.6 


±0.6 






±0.6 


±0.6 






Pion ID 










±0.2 


±0.2 


±0.2 


±0.2 


Trkeff. {K,7t) 


±1.4 


±1.4 


±2.8 


±2.8 


±2.8 


±2.8 


±4.2 


±4.2 


Ks eff. 






±0.9 


±0.9 






±0.9 


±0.9 


BB counting 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


Fisher 


±0.3 


±0.5 


±0.6 


±1.1 


±0.6 


±1.0 


±1.1 


±2.2 


BB likelihood 


±0.6 


±0.6 


±0.9 


±0.9 


±0.9 


±0.9 


±1.7 


±2.1 


Model dep. 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


±7.0 


±4.0 


±7.0 


MC statistics 


±0.4 


±0.5 


±0.4 


±0.5 


±0.5 


±0.6 


±0.5 


±0.7 


Total 


±5.2 


±5.4 


±6.3 


±6.4 


±6.2 


±8.6 


±7.5 


±9.8 



Systematic uncertainties on the signal yields obtained from the maximum-likelihood fit arise 
from three sources: uncertainties in the parameters describing the signal shapes, uncertainties in the 
combinatorial background shape, and uncertainties in the peaking backgrounds. The uncertainties 
in the means and widths of the signal shapes are obtained by comparing data and simulated data 
in charmonium control samples. For modes with electrons, we also vary the fraction of signal 
events in the tail of the AE distribution. To evaluate the uncertainty due to the background 
shape, we reevaluate the fit yields with three different parameterizations: (1) an exponential shape 
for AE, (2) a quadratic shape for AE, and (3) an uies ARGUS slope parameter C [27], which is 
linearly correlated with AE. The total systematic uncertainty in the fitted signal yield introduces 
a systematic uncertainty A;^^^ in the measured branching fraction. 

As cross checks, we also test our fit method by measuring the branching fractions and Acp 
of the J/V' /C(*)and ipi2S) K^*) final states using the vetoed charmonium events. The measured 
branching fractions are in good agreement with the 2004 world average [28] and the recent BABAR 
measurement [29]. The direct CP asymmetries Acp are all consistent with zero. We also analyze 
K^*^e^ samples and obtain signal yields consistent with zero. 

6 Results 

The results for the fits to the individual decay modes are shown in Table |31 Branching fraction 
uncertainties are predominantly statistical, with total systematic uncertainties of about 10% in 
each decay mode. 

To combine the results from the individual modes into the total B Kl^i^ and B K*i^£~ 
branching fractions, we perform a maximum-likelihood fit where the event yields in all of the modes, 
after being corrected for selection efficiency and K^*^ branching fractions, are constrained to the 
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Table 3: Results from fits to the individual decay modes. The columns from left are: decay 

mode, fitted signal yield, signal efficiency, systematic error on the selection efficiency, systematic 
error from the fit, the resulting branching fraction (with statistical and systematic errors), and the 
significance of the signal (including systematic errors). 
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K*+e+e- 




15.4 
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±0.60 


77+0-87 
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±0.60 


1.0 


B+- 




A 7+4-6 
' -3.4 


9.0 


±0.10 


±0.13 




±0.16 


1.6 



same value. In this fit we constrain the production rates of charged and neutral B meson pairs in 
the T{AS) decay to be the same. We also constrain the total width ratio T{B^) /T{B~^) to the world 
average B meson lifetime ratio T+/ro = 1.086 ± 0.017 [28]; all branching fractions from combined 
fits are expressed in terms of the B^ total width. In B ^ K*i~^£~ we perform the fit with the pole 
region included, adding the constraint: 

r{B K*n+fi-)/r{B K*e+e-) = 0.752. 

As described in Section 1, this originates from the enhanced contribution in B ^ K*e^e~ from the 
photon penguin amplitude near = 0. The branching fraction for this combined fit is expressed 
in terms of the B^ — > K*^fi^fi^ channel. We also perform combined fits to the electron and muon 
channels separately. Table 0] summarizes the results for the combined branching fractions. The 
combined significance of the signal, including statistical and systematic uncertainties, is 6.6c7 and 
5.7a for the B — > Ki~^£~ and B K*i'^i~ modes, respectively. 

The combined fits to i? — > Kl'^d." and B K*£~^i~ are shown in Figures |1] and |SJ They 
correspond to the branching fraction measurements of 

B{B Ktr ) = (0.341°;°^ ± 0.03) x 10"^ 

B{B K*tr) = (0.78i[5i? ± 0.12) X 10"^ 

where the first error is statistical and the second is systematic. The satellite peak in the IS.E 
distribution at —0.15 GeV for the B K£~^i~ fit arises from the feed-down component of the fit. 
Examination of events in this region confirms that the addition of a charged or neutral pion results 
in candidates consistent with B K*i'^i~ signal. The effect of such events on the B — > Ki~^i~ 
signal yield has been studied with fits to simulated samples, and the associated bias to the signal 
yield is negligible. 

For the combined modes we measure the direct CP asymmetries 

Acp{B+ K+i+i-) = 0.08 ± 0.22 ± 0.11 
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Table 4: Results from fits to the combined K^*^£'^£^ decay modes. The columns from left are: 
decay mode, fitted signal yield, systematic error on the selection efficiency, systematic error on 
the branching fraction introduced by the systematic error on the fitted signal yield, the result- 
ing branching fraction (with statistical and systematic errors), and the significance of the signal 
(including systematic errors). 

ABcs Afifit 13 Significance 
Mode Signal yield (lO'^) (ip-S) (ip-e) (g) 

B Ke+e~ 27.9+^;^ ±0.02 ±0.01 0.331^;^^ ± 0.02 5^3 

B Kfi+fi- 17.ll5;3 ±0.02 ±0.03 0.35i[3:i? ± 0.03 3.8 

B+ K+£+i- 36.711^ ±0.02 ±0.02 0.38t[]:[!^ ± 0.03 6.2 

^0 ^ K^i+i- S.ltii ±0.02 ±0.02 0.29t[]:i^ ± 0.03 2.8 

B Ki+i- 45.0t|^ ±0.02 ±0.02 0.34t[j:g^ ± 0.03 6.6 

B^K*e+e- ^G.ltlH ±0.06 ±0.13 0.97tJ]|? ± 0.15 4.5 

B^K*fi+fi- 20.7i^:^ ±0-08 ±0.11 0.90tJ]:i^^ ± 0.13 3.5 

B^^K*^i+r 45.2tjjj ±0.06 ±0.09 O.SltJj'ig ± 0.10 5.4 

B+ ^ K*+i+£- n^tli ±0.06 ±0.21 0.74^0-43 ± 0.22 1.5 

B^K*e+i- 56.8ti|^ ±0.05 ±0.10 0.78ti]:}^ ± 0.12 5.7 



AcpiB K*rr) = -0.03 ± 0.23 ± 0.12 

where the systematic uncertainty is dominated by the unknown asymmetry in the peaking back- 
grounds. 

Table 0] also contains the results from independent fits to the muon and electron channels, with 
no constraint enforced on the ratio of the two. From these fits we find the ratio of muon to electron 
branching fractions over the full range of to be 

Rk = 1.06 ± 0.48 ± 0.05 

Rk* = 0.93 ± 0.46 ± 0.06 

where these are expected in the Standard Model to be 1.00 and 0.75, respectively, with small 
theoretical uncertainties. 

We also perform the fit to the B — > K*£'^i^ channels with the pole region (q^ < O.lGeV^/c^) 
excluded, which modifies the Standard Model constraint on the ratio of branching fractions from 
0.752 to 1. With the pole region removed, we obtain 

B{B K*e-^e-,q^ > 0.1 GeV7c^) = (0.65lg|f ± 0.12) x 10^^ 
B{B KV+/i-,g2 > O.lGeVVc^) = (0.89+[5|^ ± 0.13) x 10"^ 
B{B K*tr,q^ > 0.1 GeVVc^) = {OMtoH ± 0.12) x 10^^ 

(1) 

From the fits to the B — > K*i'^l~ mode above the pole region, we find the ratio of muon to electron 
branching fractions 

RK*{q'^ > O.lGeVVc^) = 1.37i°:^^ ±0.11, 
which is expected to be 1.00 in the Standard Model. 
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The measured B — > K*l^i branching fraction is consistent with the previously pubhshed 
BABAR result [14] measured with 113fb-\ B{B K*£+e-) = {O.SStlH ± 0.10) x lO'^. The 
B — > K£'^£~ branching fraction is somewhat lower than the previous published BABAR result of 
B{B Ki+i-) = (0.65to;}3±0.04) x 10"^ [14]. Including correlations between events selected, the 
significance of this difference is equivalent to 2.2a. The ratios -R^(») of muon to electron branching 
fractions are also consistent with the previously published values. 

As a cross check, we have also examined the mu distribution of candidate events in the signal 
region. Of particular interest would be any evidence for a large excess in the mu spectrum near 
the lower boundaries of the veto regions, which could indicate J/^l^ or ip{2S) events escaping the 
veto. The m« spectrum, shown in Figure El exhibits no evidence for such an enhancement. The 
data points cover the full allowed region in mu, including the pole region in B ^ K*£^i^ . 




5.2 5.22 5.24 5.26 5.28 -0.2 -0.1 0.1 0.2 

(GeV/c^ ) A E (GeV) 

Figure 4: Distributions of the fit variables in K£~^£~ data (points), compared with projections of 
the combined fit (curves): (left) mss distribution after requiring —0.11 < AE < 0.05 GeV and 
(right) AE distribution after requiring |mES ~ it^b\ < 6.6 MeV/c^. The solid curve is the sum of 
all fit components, including signal; the dashed curve is the sum of all background components. 
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m^^ (GeV/ ) {Gey/c' ) 

Figure 5: Distributions of the fit variables in K*i^i'' data (points), compared with projections 
of the combined fit (curves): (left) mES after requiring —0.11 < AE < 0.05 GeV and 0.817 < 
it^Ktt < 0.967 GeV/c^, (upper right) AE after requiring |mEs — "t-bI < 6.6 MeV/c^, 0.817 < 
iTT-Kn < 0.967 GeV/c^, and (lower right) tjik-k after requiring |mEs — ^b\ < 6.6 MeV/c^ and 
—0.11 < AE < 0.05 GeV. The solid curve is the sum of all fit components, including signal; the 
dashed curve is the sum of all background components. 
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Figure 6: Distribution of mu for candidates in the signal region for B Ki~^i~ (left) and S — 
K*£~^£~ (right). The points are data. The dark gray bars are the total expectation from the sum 
of simulated signal and backgrounds, where the width of the bars reflects the uncertainty due to 
simulation statistics. The white histogram shows the signal subset of the total. 



Figure [7| summarizes the experimental measurements (points) and their theoretical predictions 
(boxes). The measurements are in general agreement with the range of rates predicted by the form 
factor calculations in Ref. [3]. The measured K^i^i^ branching fraction is significantly 

lower than the range estimated in [13]. 

7 Summary 

We have measured the branching fractions and direct CP asymmetries Aqp of the rare FCNC 
decays B K£^i~ and B K*i'^i^ . We find the (lepton-flavor-averaged, i?-charge-averaged) 
branching fractions 

B{B Ktr) = (0.34 ± 0.07 ± 0.03) x 10"*^ 

B{B K*tr) = (0.78l[!:J^ ± 0.12) X 10-^ 

consistent with the Standard Model predictions for these modes. We find Acp{B^ K^i^£~) 
and Acp{B — > K*i~^i~) consistent with zero, to a precision of 25%. We have also measured the 
ratios of the branching fractions of muon pairs to that of electron pairs; these are also consistent 
with the Standard Model to a precision of 50%. All of the measurements are statistically limited. 
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